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Recent experiments on strongly correlated bilayer quantum Hall systems strongly suggest that, 
contrary to the usual assumption, the electron spin degree of freedom is not completely frozen 
either in the quantum Hall or in the compressibles states that occur at filling factor v = 1. These 
experiments imply that the quasiparticles at v — 1 could have both spin and pseudospin textures 
i.e. they could be CP 3 skyrmions. Using a microscopic unrestricted Hartree-Fock approximation, 
CO . we compute the energy of several crystal states with spin, pseudospin and mixed spin-pseudospin 

textures around v = 1 as a function of interlayer separation d for different values of tunneling 
(Asas) , Zeeman (Az), and bias (Aj,) energies. We show that in some range of these parameters, 
crystal states involving a certain amount of spin depolarization have lower energy than the fully spin 
polarized crystals. We study this depolarization dependence on d, Asas , A-z and Aj, and discuss 
how it can lead to the fast NMR relaxation rate observed experimentally. 
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PACS numbers: 73.21.-b,73.40.-c, 73.20.Qt 



C3 . I. INTRODUCTION 



It is now well established that the two-dimensional electron gas (2DEG) in a double-quantum-wcll system (DQWS) 
at filling factor v = 1 has a broken symmetry ground state that can be described as either an easy-plane pseudospin 
q ■ ferromagnet or as an excitonic superfluid. The ferromagnetic state (in the pseudospin language) has finite interlayer 
coherence even in the absence of tunneling if the interlayer separation d is lower than a critical layer separation 
1— — '■ d c w 1.2£ where t = yJhc/eB is the magnetic length. This is an incompressible state supporting a quantum Hall effect 
I ' (QHE). The phase diagram and physical properties of this state have been extensively studied over the past fifteen 
, years (for a review, see Refs. [l|@). 
■^j- ■ In most studies of the bilayer coherent states at or near filling factor v = 1, it is generally assumed that, due to the 
strong magnetic field, the ground state is fully spin polarized. The spin degrees of freedom could thus be left out of 
the analysis. Recent experiments, however, cast some doubt on the validity of this assumption. These experiments 
include the measurement of the dependence of the activation energy of the bilayer quantum Hall state at v = 1 on 
an in- plane field, and the measurement*^ of the nuclear magnetic relaxation (NMR) time T\ near v = 1. 

The monolayer quantum Hall state at v = 1 is a spin ferromagnet. In the absence of Zeeman coupling, the lowest- 
energy charged excitation is a spin-textured topological object called a Skyrmion£i&. Measurement of the activation 
energy shows an increase of the energy gap when the magnetic field is tilted away from the z axis at v = 1. This 
is easily understood since keeping the filling factor constant means increasing the magnetic field, so the increase in 
I ■ activation energy reflects the increase in the Zeeman energy cost of the skyrmions. By contrast, measurement of 
the activation energy in the bilayer quantum Hall state at v = 1 shows a strong decrease of this energy with tilting 
until some critical angle 9 C above which the activation energy ceases to depend on 0. This anomalous behaviour of 
the activation energy was interpreted as a change in the ground state of the system at C due to a the occurence 
of a commensurate-incommensurate transition^ The lowest-energy charged excitation of the spin-polarized bilayer 
system at v — 1 is a bimeron, i.e. a skyrmion in the pseudospin field. A Hartree-Fock calculation of the behavior of 
the bimeron energy in tilted magnetic field reproduces qualitatively the features found in the experiment^ 

Sawada et al£- measured the activation energy of the 2DEG in a DQWS as a function of a parallel magnetic field 
and electrical bias between the layers. They define the imbalance parameter a — (n^ — ur) / [ul + tir) where nmi,) 
is the density in the right(left) layer of the DQWS. At a — 0, the activation energy showed the behavior expected 
for pseudospin-skyrmions (i.e., bimerons) while at a = 1, where all electrons reside in one well, the behavior was that 
expected for a spin-skyrmion. Sawada et al. found a continuous evolution from pseudospin-skyrmion to spin-skyrmion 
as the imbalance parameter was increased in various samples with tunneling energies ranging from As^s = IK to 
Asas = 33K. They concluded that the excited quasiparticles must contain both spin and pseudospin flips in order to 
explain their results. In particular, the behavior of the activation energy could not be explained by a level crossing 
between skyrmion and a bimeron excitations. Instead, the bimeron excitation, at the balance point, continuously 
transformed into a spin-skyrmion at high bias. This suggests the quasiparticles at these biases may be some object 
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that interpolates between the two types of skyrmions. Such objects have been studied in the field theoretic literature 
by Ghosh and Rajaramani^, and more recently by Ezawa and Tsitsishvilii 3 ., who dubbed these objects CP 3 skyrmions. 
In this last work, good agreement between theoretical calculations and the measurements of Sawada et al. were 
obtained. 

Another set of experiments confirming the necessity to take into account spin depolarization in the ground state of 
the 2DEG around v — 1 are those of Spielman et alA and by Kumada et al£. In these experiments, the NMR relaxation 
time Ti is measured as a function of filling factor. (In Ref. 0, this is also done as a function of an electrical bias.) 
The behavior of T\ seen in these experiments is reminescent of that measured^ in monolayer quantum Hall systems 
where the relaxation rate T± increases when v deviates from 1. A possible explanation^ involves the inclusion of 
skyrmions in the groundstate when the system is doped away from v — 1. A single skyrmion has it spin aligned with 
the Zeeman field at infinity, reversed at at the center of the skyrmions, and has nonzero XY spin components at 
intermediate distances in a vortex-like configuration. For \v — 1| > 0, the finite density of these objects is expected to 
condense into a crystal^. The quantum mean-field energy of this crystal is independent of the angle ip which defines 
the global orientation of the XY spin components. Cote et alm^ showed that this extra U(l) degree of freedom leads 
to broken symmetry and hence to a spin wave mode that remains gapless in the presence of the Zeeman field. It is 
the existence of this extra gapless spin mode in the crystal phase (and possibly in some overdamped form in a Skyrme 
liquid state) that is believed to be responsible for the rapid nuclear spin relaxation observed in the experiments. 

Our goal in this paper is to show that crystal states with some amount of spin depolarization due either to spin- 
skyrmions or CP 3 skyrmions exist around v = 1, with lower energy than crystal states with maximal spin polarization. 
We show this by comparing the energy of several crystal states in the Hartree-Fock approximation. We study the 
spin and pseudospin textures of these states as the interlayer separation, the Zeeman, tunnel coupling or the electrical 
bias are varied. Because CP 3 crystals must have a gapless spin wave mode, just as the Skyrme crystal in a monolayer 
QH system has, this crystal state could be responsible for the fast NMR relaxation rate seen in the experiments 
of Spielman et alA and Kumada et al~. In addition, these Hartree-Fock calculations provide a check on the field 
theoretic calculations that have been used on the CP 3 stat o 12 i 13 . In the latter it is assumed that the many body 
wavefunction can be represented by a local 4-spinor with constant magnitude. While Hartree-Fock introduces its own 
approximations by using a different approach we may hope to better understand the true many body state. 

Our paper is organized as follow. The model hamiltonian and the Hartree-Fock formalism needed to compute the 
CP 3 Skyrme crystal are introduced in Sec. II. In Sec. Ill, we introduce the CP 3 skyrmion and discuss its limiting 
forms: spin-skyrmion and pseudospin-skyrmion. Our numerical results are presented in Sec. IV. We discuss the 
relevance of our results to the experiments of Spielman et alX and Kumada et alA in Sec. V. 

II. MODEL HAMILTONIAN AND HARTREE-FOCK FORMALISM 

A. Hamiltonian of the 2DEG 

We consider a symmetric double-quantum-well system in a magnetic field B = Bz, and submitted to an electrical 
bias Ah = Eji — El where E R ( L ) are the subband energies in each layer (right and left) in the absence of magnetic 
field and tunneling. For the sake of limiting the number of parameters characterizing our DQWS, we make a narrow 
well approximation, i.e. we assume that the width b of the wells is small (&<<£) and treat interlayer hopping in a 
tight-binding approximation. The single-particle problem is then characterized by the separation d (from center to 
center) between the wells and the tunneling integral t = Asas/%- In the Landau gauge where the potential vector is 
taken as A = (0, Bx, 0), the Hamiltonian Hq of the non-interacting 2DEG is given by 

H Q = / 'v-- r \.,,,, r \../- (!) 

( C X,R,a c X,L,a + ^x,L, a c XM^ ■ 

X,a 

In Eq. |T|). c x . a is an operator that creates an electron with guiding center X in well j — R, L with spin index 
a = ±1. We work in the strong quantum limit where we assume that Landau level mixing is negligible and only 
the Landau level N = needs to be considered. The energies Ej iCt are defined by E^.a = Ah/2 + a.Az/2 and 
Eh.a = — Ah/2 + aAz/2 where Az = g*fisB is the Zeeman energy, g* is the effective gyromagnetic factor and /is is 
the Bohr magneton. 

We describe the various phases of the electrons in the DQWS by the set of average values I (pfj(q)) > where 
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Pi'f (l) i s an operator that we define^ by 



¥> x 



where the Landau level degeneracy is N v = S/2ir£ 2 (with S the area of the 2DEG). We explain the physical meaning 
of these operators below. 

In the Hartree-Fock approximation, the Hamiltonian of the interacting 2DEG in the DQWS is given by 



Hhf = N^ E ^P?fW 

i,ct 

a 

+ N * (5) £ E E ^ («) 0C< a (-i)> rf/ta) 



a,/3 i,J q 



where the renormalized single-particle energies E a ^ are defined by 



E ■ = E- 



2 



(2) 



In Eq. (J2J), = ^ Q and = ^ Q z^r iC[ where v i a is the filling factor for state (i, a) and v = Y]. a Vj a is the 
total filling factor of the 2DEG. The Hartree and Fock intrawell and interwell interactions are defined by 



m,i (q) 


= H(q) 


H.-Lj (q) 


= H(q) 


Xi,i (q) 


= X{q) 


x m (<?) 


= X(q) 



1 -V* 







-j-oc 

dye-y 2 / 2 Jo (qty) , 

+ 00 

dye- y2/2 e- dy/l J (q£y) . 
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B. 



Calculation of the order parameters j ^Pij^Cl)^} °f *he crystal phs 



To simplify our notation, we now define the four states 1, 2,3,4= (R, +) , (R, — ), (L, +), (L, — ) and write the order 
parameters (.Pi t 'f( f i)) as (pj,j(q)). From now on, the indices i,j will run from 1 to 4. The average values (/0i,j(q)) 
are obtained by computing the single-particle Green's function 

dj (X,X',t) = - (tc^ x (t) c\ x , (0)) , 

whose Fourier transform we define as 

Gi,i (q,r) = j r Y. ^ MX+X,) tx,x'- qy izG tJ (X,X',r) , 

so that Gi.j (q,r = O - ) = (pjj (q)). In a homogeneous phase, only {^pf/* (l = ^))} are nonzero while, in a crystal, 
(q)) 7^ only if q G {G} where {G} is the set of reciprocal lattice vectors of the crystal. 
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In our numerical calculation, we consider a finite number TV of reciprocal latttice vectors (Gi, G2, Gat). Defining 
the column vectors 

/ Gij (Gi, w„) \ 

Gij- (G2, 0J n ) 



\ (Gat, w n ) / 



where w„ is a Matsubara frequency, and the vectors 



f 1 ^ 
















,0 = 




W 




W 



along with the 47V x 4 matrices 



G 



f Gi,i G_i,2 G_i,3 Gi^ \ 

G2.I G_2,3 G_2A 

G3.1 G3.2 G_3 : 3 G_3,4 

\ G4.1 G4.2 G4.4 J 



and 



C = h 



I B _0 \ 

B 0_ 

B _0 

V ob/ 



we find that the Hartree-Fock equation of motion for the single-particle Green's function matrix G can be written in 
a matrix form as 



(iftwn + p)G- AG = C, 



(3) 



where A is the 47V x 47V hcrmitian matrix 



A17 
-X (pit) 7 
i^G,G^ — X (pi, 3) 7 
-*<PM>7 



with 7 -e-< GxG ')^ 2 /2 7 



-X(p 2 ,l}^ ~t^G,G' ~ X (p3,l) 7 -^(P4,l)7 \ 



^A 2 7 

(P2,3) 7 



(P3,2) 7 
A37 



-^G,G' - -X" (P2,4) 7 -X{p 3A )j 



and 



-X (p 4 ,3>7 
A47 



T 4 = Ti (G - G') = £ (G - G) ( aj (G - G')> . 



(4) 



In the definition of A, the symbols Aj, X and AT stand for A 4 (G — G') , X (G — G') , TV (G — G') and the quantity 
(p i , j )for(p i , j (G-G / )). 
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The (pi : j (G))'s are found by solving numerically the self-consistent equation of motion given by Eq. This 
equation has many solutions representing the local minima of the Hartree-Fock energy given by 



E K^ +2 (0)) + (p l+ 2,(0)>] 



^ (5) E E ^ ( G ) (- G )) <^ («)) 



i,j G#0 

1 /e 2 



2zy V Kt 



There is no guarantee that a solution is an absolute minimum of Ehf- Instead, we compare a finite number of 
likely solutions and choose the one that minimizes the energy. The numerical scheme to solve Eq.Q is described in 
more detail in Rcf. 

By definition, 

E fa.< (°)> = 

i 

This equation fixes the chemical potential fx in Eq.|j3J). 



C. Spin and pseudospin fields in the crystal phases 

In the Landau gauge and with the Hilbert space restricted to the first Landau level only, an electronic state in a 
single quantum well system (SQWS) is specified by the two-component spinor cx where 



cx 



c x,+ 
cx,- 



Similarly, an electronic state in a spin-polarized DQWS can be described by mapping this two-level system into a spin 
1/2 system by using the pseudo-spinor cx where 



, c x,R 
c x = „ 

cx,l 



For a four-level system (with states j = 1,2,3,4= (R, +) , (R, — ), (L, +), (L, — ) as given above), an electronic state 
is specified by the four-component spinor 



cx 




The operators Pi,j(q) that we introduced previously can be mapped into the density operator, p(q), the spin and 
pseudospin densities operators S a (q) and P a (q) (with a = x,y, z) and the 9 operators R a ,b (q) using the SU(4) 
algebra (as defined in Ref. IT^) 

= ^E e ^ +lM//2 *W-**> (5) 
^ v x 



x 



,-iq x X+iq x q v i 2 /2 t 



^X'a 



Cx- 



■1v 



(G) 
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P*(q) = ±^2e-^ x +^ 2 / 2 c^ x rr n c x _ qyP , 



(7) 



x 



(8) 



where the 4x4 matrices r| pm and T^ pm are defined by 



with a a a Pauli matrix, and by 



r spin 



<Ja 
<J a I ' 



pptn 



? 
















o 













where / is the 2x2 unit matrix. 

From Eqs. (JS])-©, it is easy to show that the electronic densities in the right and left wells are given by 

(PR (q)> = (Pi,i (q)> + (P2,2 (q)) , 

(pl (q)> = (p3,3 (q)> + (PiA (q)> , 

with the total electronic density given by (p (q)) = (pn (q)) + {pl (q)) • 
The spin densities in the right and left wells are given by 

(S x , R (q)) = R(pi, 2 (q)), 
(5 VlK (q)) = 3(pi. 2 (q)>, 

(S^(q)) = 5R(p 3 , 4 (q)>, 
(q)> = 3 (p 3 ,4 (q)) , 

and by 

(^(q)> = 5 [<Pi,i (q)> - <P2,2 (q))] , 

(S,,M) = i[<P3,3(q))-<P4, 4 (q))]. 

Finally, the pseudospin densities for the up (+) and down (-) spin components are given by 

(P x ,+ (q)) = »<Pi,s(cO>, 
(P„,+ (q)) = 9f<Pi, 3 (q)>, 

(P x ,-(q)> = 5R<P2, 4 (q)), 
(P„,-(q)) - S(p 2 , 4 (q)>, 

and by 



(P,,+ (q)) = 5 [<Pi,i(q)>-<P3,3(q))], 
(P,,-(q)) = i[(p 2 , 2 (q))-(p 4 , 4 (q))], 



with P z (q) = P z ,+ (q) + P Z) - (q) the total pseudospin density. 
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Note that by definition {p id (q)) = (p jti (-q))* . Also, (S (-q)> = (S(q))* and (P (-q)) = (P(q)>*. The four 
order parameters that are not related to the electron, spin, or pseudospin densities are (pi^ (q)) , (P2.3 (o.)) an d their 
complex conjugates. These densities involve average values of operators that flip both the spin and the pseudospin. 

The Hartree-Fock energy per electron can now be written as 



Erf 
N 
1 



^ (P z (0)> - ^ (S z (0)) - (P x (0)) 



j^J2 Ti ( G ) 1 (p ( G ))i 2 + J E j m ( G ) 1 ^ ( G )) 

G G 

-JEE A(G)|(S a (G))| 2 

G a=R,L 

;E E ^(G)[|(F^(G))| 2 + |(P a , Q (G))| 2 

G a=+,- 

-I^l(G) [|(p M (G))| 2 + |( P2 , 3 (G))| 2 



(9) 



In Eq. @, we have defined the interactions 



T X (G) = H(G)+H(G)--X(G), 
J Ztl (G) = H(G)-H(G)-±X(G). 
Because of the neutrality of the total system comprising the electrons and the positive donors, we have 

Tr(0) = -~X(0), 
^,i(0) = j-\x{Q). 

In the pseudospin langage, a bias acts as a pseudo-magnetic field that couples to the z component of the total 
pseudospin while the tunneling acts as a pseudo-magnetic field that couples to the x component of the total pseudospin. 
The positive sign in front of the first term on the r.h.s of Eq. JSJ is due to our particular choice of mapping (R — > + 
and L — ► — ) for the pseudospin states. A positive bias forces the z component of the pseudospin down; i.e. pushes the 
electronic charge in the left well. At zero bias, there is equal population of electrons in both wells and (P z (0)) — in 
order to minimize the capacitive energy. 



D. The coherent liquid state at v — 1 



At v = 1, the 2DEG can have spontaneous interlayer coherence. For nonzero Zeeman coupling, the ground state 
is well described (but this description is only exact at d = 0) by a state where all electrons are in the symmetric 
combination of both wells and all spins are polarized. The order parameters are then 

(Pl,l (0)) = (p 3 ,3 (0)> = \, 

(Pi,3(0)) - \e l \ 

irrespective of the value of d (for t 7^ 0, 8 = where 9 is the angle between the pseudospin vector and the x axis). 
In the absence of tunneling, the coherent liquid phase supports a gapless pseudospin wave excitation^ 8 that disperses 
linearly with q (for d ^ 0) at small wavectors and becomes soft at an interlayer separation d c /£ » 1.2. This critical 
separation is increased by a finite tunneling. Above this critical interlayer separation, the interwell coherence is lost. 
The system is then believed to be formed of two composite fermions liquids with filling factor vr = vl = 1/2. This 
state is not captured by the HFA which instead predicts a transition to a charge-density-wave state. 
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E. Influence of a bias 

With a bias, the symmetric and antisymmetric states of the non-interacting 2DEG are replaced by the bonding 
and anti-bonding states defined by 

where a = . Ab = is the unbalance parameter. At v = 1, the ground state has all electrons in the IB) state (i.e. 

V A l+ A SAS 

the symmetric state in this case) with up spin if the Zeeman coupling is non zero. 

When interactions are included, we can still easily solve Eq. J3J) in the presence of an electric bias, at v = 1, 
assuming that the Zeeman term is non zero so that the 2DEG remains spin polarized. We find 

(Pis (0)) = v R = ~ (1 - a) , 

(P3,3(0)) = I/ £ = l(l+(7), 

(Pi,a(0)) = a = Vv^. 

When t = 0, the energy of the 2DEG is again invariant with respect to 9 while for t ^ 0, the energy is minimized 
when 9 = 0. The bias acts as a pseudomagnetic field that forces the pscudospin up or down from the xy plane. The 
interlayer coherence is maintained but diminished. In the HFA, there is a critical interlayer separation d^ F (A&) where 
interlayer coherence is lost and all the charge is transferred into one well. Notice that the pseudospin mode remains 
gapless under bias although it now becomes soft at a critical interlayer separation d^ RPA (A&) < d^ F (A&) that 
depends on bias. This situation is represented in Fig. 1 in the case of zero tunneling (essentially the same calculation 
can be done at nonzero tunneling. As expected, the critical interlayer separation increases with As^s)- For a 2DEG 
initially in the incoherent state at d > d® RPA (A&), it is possible to get a coherent state (and so a quantum Hall 
effect) by increasing the bias. This transition has been studied in detail both theoretically^ and experimentally^. 
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FIG. 1: Phase diagram of the 2DEG in a DQWS at v = 1. The top critical line (triangles) is obtained by the instability of 
the pseudospin wave mode in the GRPA. As d/£ goes through this line, the interlayer coherent state (ICS) loses its coherence 
to become two interlayer incoherent states (IIS) with Vril) = 1/2- When bias is increased at fixed d/£, the ICS looses its 
coherence and all the charge is transferred into a single quantum well (SQW liquid). This line (squares) is obtained in the 
HFA. 



III. TOPOLOGICAL EXCITATIONS IN THE DQWS 
A. Spin and pseudospin skyrmions 

In a monolayer system with only spin degrees of freedom, Eq. becomes 
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-L^t 2 (g)|(p(g))| 

G 

-i^A(G)|(S(G))| 2 , 

where the interaction 



^ = -^(0)) (10) 
TV v 

1 



V 
G 



T 2 (G) = 2H (G) — X (G) . 

In the absence of Coulomb electrostatic energy and Zeeman coupling and in the gradient approximation where 
(S (r)) is assumed to vary smoothly over the magnetic length i, the energy in Eq. (|10f> can be reduced to that of the 
0(3) nonlinear sigma models The 0(3) nonlinear sigma model in a planar geometry possesses topological solitons or 
skyrmions. Using the complex function w (z) which represents the stereographic projection of the unit sphere of spin 
textures s, a skyrmion of Pontryagin index Q = 1 (which corresponds to the addition of one electron to the 2DEG) 
can be written^ as 

l-s z (z) A 

where z = x + iy. Eq. I|ll|) describes a skyrmion of size |A| centered at position b = x^ + iyb in the x — y plane. The 
spin components are 



2A* (z - b) 
Sx(Z) A) ~ |A| 2 + |z-&f 

\z-b\ 2 + \X\ 2 

At infinity, the spins point upward while at the center z = b of the skyrmion, they point downward. The phase tp in 
A = |A| e lip fixes the global orientation of the spins forming the skyrmion. 

In quantum Hall systems, skyrmion-antiskyrmion excitations have been shown to have lower energy than the corre- 
sponding maximally spin-polarized Hartree-Fock electron-hole excitation if the Zeeman coupling is not too strong£22 
Around filling factor v = 1, a finite density of skyrmions (y > 1) or antiskyrmions (y < 1) are included in the 
groundstate (with filling fraction v s = \ v — 1|). At T = K, these quasiparticles condense into a Skyrme crystal. 

The equation of motion method derived in Ref. ^3 was used some years ago to study the phase diagram of the 
Skyrme crystal in the Az — v s space. Numerical result o 1 ^' 1 ^ show that, in a large portion of the phase space, skyrmions 
crystallize into a square lattice structure with two skyrmions of opposite phases (ip — and ip = it) per unit cell. 
This arrangement was termed SLA (for square lattice antiferromagnetic). More recent calculations^ show that, as 
the Zeeman coupling is increased from small values at fixed quasiparticle filling v s , the crystal structure changes from 
a checkerboard 2 ^ lattice of merons near Az — (i.e. a SLA skyrmion lattice where each skyrmion splits into two 
merons of charge e/2 with opposite vorticities and where all merons are equally spaced) to a lattice of biskyrmions 2 ^ at 
small Az, a SLA skyrmion lattice at moderate Az and finally into a TL120 lattice of skyrmions (a triangular lattice 
with three skyrmions with phases <p = 0, 27r/3, 47r/3 per unit cell to avoid the frustration created by the preferred 
antiferromagnetic ordering of the skyrmions^) at higher values of Az- The spin texture is gradually lost as the 
Zeeman coupling further increases and we finally have a Wigner crystal of maximally polarized quasiparticles with no 
spin texture. 

For a spin-polarized 2DEG in a DQWS, the energy functional of Eq. becomes 

Ehf £ (P z (0)) - ^ (P. (0)) 



N 



^£T 2 (G)|(p(G))| 2 (12) 

G 

IJ2 J ^ (G)|(P,(G))| 2 

G 

TIE E ^(G)|(P^(G))| 2 , 



V 

G a=+, 
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with 

J z , 2 (G) = H (G) - H (G) - X (G) . 

In the absence of bias, tunneling and Coulomb electrostatic energies and in the gradient approximation, the energy 
in Eq. can be reduced to that of the anisotropic nonlinear sigma model with a unit pseudospin field p (r)i2i The 
topological excitations of this model are bimerons and merons (or pseudospin-skyrmions and pseudospin- merons). A 
bimcron with topological charge Q — 1 has its pseudospin field given bj^ 



p x (z) - ip y (z) 
Pz (z) 



2 - b) (z* + b*) 
\z-b\ 2 + \z + b\ 2 '' 
\z - b\ 2 - \z + b\ 2 

1 2 . i . ,i2' 



\z-b\ z + \z + b\< 

where +b are the positions of the two merons forming the bimeron. Alternatively, we can write 

p x (z) - ip y (z) 



w (z) 



1 ~Pz (z) 
z - z L 

z - z R 
z-b 



z + b 



where ZutL) are the positions of the merons in the right and left wells. The angle ip gives the global orientation of the 
pseudospin vectors with respect to the x axis at infinity. When z = b(—b), we are at the center of the meron on the 
left (right) well and there the pseudospin p z = — 1(+1). 

Bimerons and merons have been studied extensively in the context of the QHE.iiSSiSi Again, bimeron-antibimeron 
excitations have been shown to be the relevant excitations near v = 1. Although we have not performed an exhaustive 
calculation of the phase diagram of bimeron crystals, our numerical results^ show that at finite tunneling, an SLA 
(or rectangular antiferromagnetic) configuration of bimerons is stable up to an interlayer separation djl w 1. At 
very small tunneling, the bimeron lattice becomes a lattice of merons with again the checkerboard configuration. In 
comparison with Wigner crystal in bilayer systems where a one-component Wigner crystal can only be stabilized at 
small interlayer distances of order d/£ ~ 0.1, the interlayer coherence in a bimeron or meron lattice persists to much 
larger d/L 



B. CP 3 skyrmion 

When both spin and pseudospin are active degrees of freedom, electronic states must be described by a four- 
component spinor. As explained in Ref. IT2I this spinor is a CP 3 spinor since the DQWS has a U(l) gauge invariance 
(all four components of the spinor must be transformed by the same phase for the DQWS's energy to remain invariant). 
Strictly speaking, a texture of a CP 3 spinor can lie wholly in the spin degrees of freedom, or wholly in the pseudospin 
degrees of freedom. We only have a guarantee that the topological charge associated with the texture integrates to an 
integer. In this paper we use the phrase "CP 3 Skyrmion" to refer to textures in which both the spin and pseudospin 
degrees of freedom are appreciably varying. This is sometimes referred to as an interwined texture. 

To start the iteration process needed to solve Eq. Q, we must provide an approximate solution for the crystal 
of CP 3 skyrmions. From our discussion above, we expect that an SLA configuration of skyrmions could be a likely 
solution. Following RajaramaniS, we write the four-component spinor 



\ A 2 

were z = x + iy and the normalisation factor is 



^|A 1 | 2 + |A 2 | 2 + 2|z| 2 + |6| 2 
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This state contains a skyrmion of size |Ai| at position b in the right well, a skyrmion of size | A2 1 at position —b in the 
left well and a bimeron in the spin up component of the pseudospin centered at z = 0. The CP 3 static energy of this 
skyrmion is given by 

e cp3 = 2 Ps Y / Yl J dr (^-*« ( r ))* ( r )) . ( 14 ) 

a=lj=x,y 

where Dj = dj — iKj with K a gauge field defined by Kj — (r)%*ii (r) and p s a "spin-pseudospin" 

stiffness. The field $ (r) must satisfy the constraint ^ a ( r )| 2 = 1- Eq. can ^ e obtained from our Hartree- 
Fock energy by taking the limit Ag^s = Az = d = 0. The solution of Eq. 113|) is a skyrmion with topological charge 
Q = J dr5Q (r) = 1 . The CP 3 topological charge density is defined by 

SQ (r) = ~i £y ' (a * q (r)) * {D ^ a (r)) ■ (15) 

The order parameters for this single quasiparticle state are given, in real space, by (pij (r)) = (r) \tj (r)^ . 

Fourier-transforming this expression, we can easily write the (pij (q)) s for this state. To write an approximate 
solution for a crystal of these quasiparticles, we consider the change 5 (pij (r)) in the ground state (at v — 1) when a 
skyrmion is added to the system 

s (pi,j (r)) = (pi,j (r)) - (pij (r)) GS , 

where (pij (r)) Gfi . describe the ground state (at v = 1) which has (at zero bias) all electrons in the bonding (or 
symmetric) state with up spins (see Sec. II). In principle, the fields S {pij (r)) are zero when we are far away from 
the position of the skyrmion so that the crystal state can be written approximately as 

(p id (r)) = ( Piij (r)) GS + E 6 (P% (' " R " Ca)) , 

R a=l,2 

where c a is the position vector of the two skyrmions in the unit cell, R is a lattice vector, and a is the index of the 
phase of each of the two skyrmions. The order parameters for the crystal state are then given by 

{Pij (G ^0)>~ J2 e- iG ' c ^(pg 3 (G)} . (16) 

a=l,2 

In general, it is sufficient to give the (pi.j (G 7^ 0))' s given by Eq. (|16fl on the first or first two shells of reciprocal 
lattice vectors in order for the program to converge to the CP 3 skyrmion crystal. The SLA configuration is obtained 
by choosing R4 = nax and R2 = may for the lattice vectors (n, m = 0, ±1, ±2, ...), Ci = — f y and C2 = +f y for the 
positions of the two skyrmions. For the first skyrmion, we take b = — | and Ai = A2 = 1. Eor the second skyrmion, 
we take, b' = | and \[ = \' 2 = —1 in order to rotate the global phase by tt. 

IV. NUMERICAL RESULTS ON CP 3 CRYSTALS 

A. Crystal states considered 

The HFA does not contain the correlations necessary to describe the ground state at djl > 1.2 where interwell 
coherence is lost. For this reason, we limit our numerical calculations of crystal states to interlayer separations 
< d/t < 1.2. In the monolayer and polarized bilayer limits, we found that a square lattice with two skyrmions of 
opposite phases per unit cell is the ground state in a wide region of parameter space. We thus choose to consider the 
following states in our analysis: 

• CP 3 : a square lattice with two spin-pseudospin skyrmions of opposite phases per unit cell as described at the 
end of Sec. III. This state is represented in Fig. 2 in the case of small tunneling where each skyrmion is broken 
into two merons of opposite vorticities. 
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• SPB: a spin-polarized square lattice with two bimerons of opposite phases per unit cell. The spinor of Eq. 

/ * - h \ 

(|13fl is replaced by (r) = A z + b I ' ^ smau t unne hng, the bimerons split into a pair of two merons 

V o 7 

with opposite vorticities. This spin-polarized bimeron (or merons) lattice should be the ground state when the 
Zeeman energy is of the order or bigger than the tunneling energy. 

• SS: this is a symmetric skyrmion state which is a pseudospin-polarized square lattice state with two symmetric- 
band spin-skyrmions of opposite phases per unit cell. By "symmetric-band", we mean that the SU(2) spinor 

I * 

for the first electron would be given by \& (r) — A' 1 




V 



in the symmetric- antisymmetric basis (S+; 



S-;AS+;AS-) basis or by * (r) = A 




in the (R+; R—;L+;L—) basis. We expect this phase to 



be the ground state state when tunneling energy dominates the Zeeman energy. 

• High Tunneling CP 3 (HCP 3 ): a square lattice with two spin-pseudospin skyrmion of opposite phase per unit cell. 
The difference between this state and the CP 3 state above is that here the spin texture exists in the symmetric 
and antisymmetric bands while in the latter it exists separately in each quantum well. This state is the ground 
state only when the tunneling energy is higher than the Zeeman energy and only for filling factor v > 1. We 
note that the HCP 3 state is an intermediate state between SS and SPB states. The textures in the HCP 3 state 
splits into two vortices with charge e/2 by reducing the Zeeman gap. 

When the tunneling or Zeeman couplings increase, the size of the pseudospin or spin skyrmions decreases. At large 
values of both these parameters, a limit is reached where the skyrmion state reaches maximal spin and pseudospin 
polarization. The resulting state may be viewed as a crystal of Hartree-Fock quasiparticles. When interlayer coherence 
is non zero, the HF quasiparticles are delocalized in both wells and form a coherent Wigner crystal (i.e. a "one- 
component "Wigner crystal) 




FIG. 2: (Color online). Spin textures in each well (bottom) and pseudospin textures in each spin component (top) of a CP 3 
crystal for v = 0.8, A SAS / (e 2 /^) = 0.0002, A z / (e 2 /n£) = 0.002, d/l = 1.0. The crystal has four merons per unit cell. In 
each unit cell, two merons with the same vorticities have opposite phases as explained in the text. The contour color indicated 
on the legends at the right side of each plot is for the z component of each field. 



With increasing bias, the CP 3 or SS solution will continuously transform into a monolayer spin skyrmion. The SS 



solution, for example can be written in the R/L basis as \& (r) = 
a is the unbalance parameter and a = in the absence of bias). 



( VT~^ z \ 

V T~~q : A i 

\fTToz 



in the presence of a bias (where 



B. Effect of interlayer separation at zero bias 



We first compute the energy of the three states just introduced as a function of the interlayer separation for 
parameters v = 0.8, As as/ (e 2 / ' nt) = 0.0002 and Az / (e 2 / 'nl) — 0.002. Figure 3(a) shows the energy differences 
Ecp 3 — Espb and Ess — Espb- At small interlayer separation, the ground state is the SPB crystal while above 
some critical interwell separation that depends on the Zeeman coupling, a CP 3 crystal state emerges and remains the 
ground state up to the largest value of djt that we consider (i.e. djt — 1.2). Figure 3(b) shows the difference in 
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energy between the CP 3 and the SPB crystals for several values of the Zeeman coupling. As expected, the interlayer 
separation at which the SPB-CP 3 transition takes place increases with increasing Zeeman coupling. (We remark that 
e 2 /n£ = 50.489\/S K, so that the difference in energy, at d/l = 1.0 and A z / (e 2 /nl) = 0.002, is of the order of 
90 mK.) Moreover, in Fig. 3(c), we see that the spin polarization per electron (S z (q = 0)) jv varies strongly with 
interlayer separation in the CP 3 crystal state in comparison with the SS state. The spin depolarization of the CP 3 
state increases with decreasing Az and reaches a maximum at about d/£ ss 1. Figure 3(d) shows the behaviour of 
the pseudospin polarization per electron in the x — y plane i.e. (P x (q = 0)) jv. The pseudospin polarization increases 
when the spin polarization decreases. The value of pseudospin polarization gives some indication of the size of the 
bimcron in the CP 3 skyrmions. When (P x (q = 0)) jv = 0.5, there are no pseudospin vortices in the phase considered. 
Note that we have chosen in our analysis a very small value of the tunneling constant. Our results of this section stay 
essentially the same if the tunneling coupling is exactly zero. 

The observation that this CP 3 state is most prominent at large d suggests that it is stabilized by the interlayer 
charging energy. The merons of the SPB state involve "tilting" of the pseudospin at their cores into one layer or the 
other, at an energy cost of order e 2 d/ 1 n£ 2 . For large enough d, it is energetically favorable to admix spin states so 
that near their cores the charge of the vortices will be balanced. An examination of the charge densities in each well 
reveals that the CP 3 lattice is indeed more locally balanced than the SPB lattice. 

C. Effect of tunneling at zero bias 

In Fig. 4, we show the behavior of the spin polarization per electron with filling factor for v < 1 for two values of 
the interlayer separation. At v = 1, (S z (q = 0)) /v = 0.5. Away from v = 1, the spin polarization decreases rapidly 
for the CP 3 crystal until it reaches a minimum at about v = 0.9. Then, as v is further decreased and the density of 
skyrmions increases, the size of the skyrmions also decreases and the Wigner crystal limit is reached where the ground 
state is again fully spin polarized. The behavior of the spin polarization we find for the CP 3 crystal is identical to 
what was found for Skyrme crystals in a single layer system^ As v is increased towards v = 1, we also find that the 
critical interlayer separation for the transition between the SPB and the CP 3 states decreases so that the CP 3 crystal 
state is stable over a larger range of interlayer separation for smaller quasiparticle filling. For the SS state that occurs 
at high tunneling, the variation of the spin polarization is less marked than in the CP 3 crystal. 

Because the spin polarization is minimal around d/£ = 1 for v = 0.8, we choose this value of the interlayer separation 
to study the effect of tunneling on the spin polarization. Figure 5 shows the difference in energy Eq P 3 — E$pb and 
Ess ~ Espb for three values of the Zeeman coupling. At small Zeeman coupling, increasing Asas causes a transition 
from the CP 3 to the SS crystal. At larger Zeeman couplings, where the ground state is the SPB crystal at zero 
tunneling, increasing As as causes first a transition to a CP 3 crystal (in a very narrow range of As as) and then 
into a SS at larger tunneling values. The CP 3 crystal exists only in narrow range of Asas and that range decreases 
with increasing Zeeman coupling so that the CP 3 state disappears at large Az- Typically, the SS phase occurs for 
A S as > A z /2. 

The change in the spin and pseudospin polarizations of the CP 3 and SS crystals with Asas is shown in Fig. 6 
for the values of the Zeeman coupling considered in Fig. 5. We see from this figure that increasing Asas increases 
the pseudospin polarization and decreases the spin polarization. At the transition from the CP 3 to the SS state 
indicated by the vertical dashed lines in Fig. 6, there is a sharp drop in the spin polarization. This sudden change in 
(S z (q = 0)) , and so in the in-plane spin polarization, should lead to abrupt changes in the NMR relaxation time. 

Our results, so far, have been for filling factor v < 1. Interestingly, we do not find a precisely analogous spin- 
pseudospin configuration for filling factors v > 1. Instead, as explained in Section IV, we have found an intermediate 
spin-pseudospin state at large values of As as and small separations ( d/t < 0.7 ), which we call the HCP 3 state. 
Fig. 7(a) shows the difference in energy of HCP 3 or SS and SPB. As we can see in this figure, by increasing Asas 
the ground state changes from SPB to HCP 3 and then to SS. Also in Fig. 7(b) we can see the spin depolarization in 
HCP 3 state as a function of Asas- The spin polarization of the HCP 3 state is more sensitive to Asas than the SS 
state. 

For v > 1 the SS state also can be the ground state. At v = 1.04, we find that the SS state is the ground state 
for djl < 0.8, As as I {e 2 / k£) > .03 and < Az/ (e 2 / nt) < 0.002. The spin polarization with interlayer separation 
in the SS state is shown in Fig. 8. The linear behaviour is typical of what is obtained for v < 1 in the SS state (see 
Fig. 3, for example). 
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FIG. 3: (Color online). CP 3 crystal at v = 0.8 and A S as = 0.0002 [e 2 /nt). (a) Energy difference between the CP 3 or SS state 
and the SPB state at A z = 0.002 (e 2 /^); (b) energy difference between the CP 3 and SPB state for different values of the 
Zeeman coupling; (c) Spin polarization and (d) pseudospin polarization per electron in the CP 3 state (full lines) and SS state 
(dashed lines) for the different Zeeman couplings indicated in (b). 
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FIG. 4: (Color online). Spin polarization in the CP 3 and SS skyrmion crystals as a function of filling factor for v < 1. Here 
A z / (e 2 /K,£) = 0.006. For the CP 3 curves, As as/ (e 2 / 'k£) = 0.0002,while for the SS curve, As as/ (e 2 / 'nl) = 0.05. 



D. Effect of a bias on the spin polarization 



To conclude this section, we look at the effect of a potential bias on the spin polarization. Intuitively we understand 
that a CP 3 skyrmion involves a "twist" in some high dimensional space that is difficult to plot in 2D. That twist will 
occur through degrees of freedom where it costs the least energy, and the texture will vary slowly in sectors where 
the system is "stiff" . If we can change the stiffness of textures along some direction of phase space we can drive the 
texture into or out of that degree of freedom. A simple analogy would be to drive an 0(3) model into an XY model 
by making excursions into the z-direction too costly. 

The behavior of this system with respect to bias illustrates this physics. We choose the parameters d/i = 1, v = 
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FIG. 5: (Color online). Energy difference between the CP 3 crystal (filled squares) or SS crystal (filled triangles) and the SPB 
state at v = 0.8 and d/£ = 1.0 as a function of tunneling for Zeeman couplings: (a) Az / (e 2 / 'k£) = 0.004;(b) Az / (e 2 / k£) = 
0.006; and (c) A z / (e 2 / 'nl) = 0.008. 



0.8, Ag as/ (e 2 /«£) = 0.0002 and Az/ (e 2 / k£\ = 0.01 so that, at the balanced point, the ground state is a spin- 
polarized meron crystal (SPB). Our numerical results, plotted in Fig. 9, show that there is a transition first into a 
CP 3 crystal and then into the SS state as the applied bias increases. The energy of the CP 3 crystal interpolates nicely 
between the SPB and SS phases as can be seen in the figure. The corresponding spin polarizations are shown in Fig. 
9(b). The bias has the effect of inducing a linear spin depolarization of the 2DEG in the CP 3 state. In effect, the 
texture inducing the deviation of charge density from v = 1 is being shifted from the pseudospin degree of freedom 
to the spin degree of freedom in a continuous way. Note that the spin polarization varies only slightly with d/t in the 
SS state. Fig. 9(c) shows the filling factor vr and vl in the CP 3 state (the exact same curves are obtained in the SS 
state). Above A5/ (e 2 / ' k£) ss 0.30, all the charge is transferred into the right layer and the spin polarization is that 
appropriate for a monolayer Skyrmion crystal with filling factor v — 0.8 and Zeeman coupling Az/ [e 2 /nltj = 0.01 
and is then independent of the interlayer separation. We expect the marked decrease in the spin polarization with 
bias to translate into an increase of the NMR relaxation rate with increasing bias. 
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FIG. 6: (Color online). Spin (a) and pseudospin (b) polarization per electron as a function of tunneling for the CP 3 crystal 
(lines with symbols) and the SS state (dashed lines) for different values of the Zeeman coupling. In (a), the spin polarization 
goes to that of the SS state at large tunneling. The vertical lines in (a) indicate the critical tunneling for the transition between 
the CP 3 and the SS states as found from Fig. 5. 



V. DISCUSSION AND CONCLUSION 



Our numerical calculations show that crystals involving spin and/or pseudospin textures are likely candidates for 
the ground state of the 2DEG in a bilayer quantum Hall system around filling factor v = 1. At small tunneling and 
for v < 1, we find intertwined spin and pseudospin textures (CP 3 crystal) with a spin polarization that is strongly 
interlayer dependent while at higher tunneling, a symmetric skyrmion state with fully polarized pseudospins or another 
type of spin-pseudospin state minimizes the energy. 

As mentioned in our Introduction, a Skyrmion crystal has an extra gapless spin mode in the crystal phase (and 
possibly in some overdamped form in a Skyrme liquid state) that is believed to be responsible for the rapid nuclear 
spin relaxation observed in the experiments This extra Goldstone mode is present both in the SS and CP 3 crystal 
states that we studied in this paper but not, in the SPB stated To make a direct comparison with the experiments 
of Spielman et al. and Kumada et al., it is necessary to compute the NMR relaxation rate. Results of such calculations 
will be presented elsewhere^ We can expect, however, that the relaxation rate will be proportional to the in-plane 
spin polarization so that the behavior of spin polarization S z should be an indication of the behavior of the relaxation 
time T]_. The rapid change in the spin polarization that we found in the CP 3 crystal state (very small Ag^g) for 
filling factor around v — 1 may explain the rapid change in the NMR relaxation rate measured in the experiment of 
Spielman et al. which was carried on at almost zero tunneling and for a Zeeman coupling which is approximately that 
indicated in Fig. 4. 

Our Hartree-Fock calculation indicates that the ground state at higher tunneling is a SS state instead of a CP 3 
crystal. In this case, the spin polarization S z varies much less rapidly with filling factor than for CP 3 crystal. Moreover, 
the spin polarization does not depend much on the interlayer separation as can be seen, for example, in Figs. 3. The 
results of Kumada et al. showing a rapid change in the NMR relaxation rate as well as a strong dependence on the 
interlayer separation would be more readily explained by a CP 3 crystal state than by the SS state that we find. This 
is also true for their measurement of the relaxation rate in the presence of an applied bias. 

Fig. 10 illustrates the in-plane component for various Zeeman couplings as a function of bias, which we believe is a 
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A SAS /(e 2 /Kl) 



FIG. 7: (Color online), (a) Energy difference between the HCP 3 crystal or SS crystals and the SPB state as a function of 
tunneling for v — 1.2, Az / (e 2 / k£) = 0.0015 and d/l = 0.1; (b) Spin and pseudospin polarization as a function of tunneling for 
the HCP 3 and SS crystals. 




d/l 



FIG. 8: (Color online). Spin polarization as a function of interlayer separation in the SS state for different values of the Zeeman 
coupling and filling factors. Here Asas/ (e 2 / k£) = 0.04 and v = 1.04. 



measure of the NMR relaxation rate, for small A^s- The evident continuous behavior is reminiscent of the Kumada 
results. We speculate that the effects of finite well width, Landau level mixing, and possibly disorder, all not included 
in our calculations, may stabilize the CP 3 state at higher tunneling than was found for our idealized system. 

In conclusion, we have studied textured quantum Hall states in bilayer systems including the spin degree of freedom, 
and have shown under appropriate circumstances that mixed spin-pseudospin textures appear as the groundstate. 
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FIG. 9: (Color online), (a) Energy difference between the CP 3 or SS states and the SPB state; (b) spin polarization per electron 
as a function of bias for v = 0.8, As as/ (e 2 /k£) — 0.0002, and Az / (e 2 / 'k£\ = 0.01; and (c) filling factor in the right and left 
wells in the CP 3 state. 
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FIG. 10: (Color online). Average in- plane spin polarization in the CP 3 crystal state as a function of applied bias for v = 
0.8, AsasI (e 2 /Ki) = 0.0002 and d/t = 1.0 
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